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ABSTRACT 

Interpretation of He n Lya absorption spectra after the epoch of He n reionization requires 
knowledge of the He n ionizing background. While past work has modelled the evolution of 
the average background, the standard cosmological radiative transfer technique assumes a 
uniform radiation field despite the discrete nature of the (rare) bright quasars that dominate 
the background. We implement a cosmological radiative transfer model that includes the most 
recent constraints on the ionizing spectra and luminosity function of quasars and the distri- 
bution of IGM absorbers. We also estimate, for the first time, the effects of fluctuations on 
the evolving continuum opacity of the IGM. Our model results in a He n ionizing background 
that evolves steeply with redshift, increasing by a factor > 3.5 from z = 3.5 to z = 2.5. This 
causes rapid evolution in the mean He n Lya optical depth - as recently observed - without 
appealing to the reionization of He n. Such behaviour could instead result from rapid evo- 
lution in the mean free path of ionizing photons as the helium in higher H i column density 
absorbers becomes fully ionized. 

Key words: cosmology: theory - intergalactic medium - diffuse radiation 
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1 INTRODUCTION 

The ionizing background is crucial for understanding many aspects 
of large-scale structure and galaxy formation at high redshifts. For 
example, unraveling the physical density structure of the Lya for- 
est (which contains most of the the intergalactic medium, or IGM, 
at z > 2) requir es knowledge of the ionizat ion state of the inter- 
vening material dRauchll 19981 ; lMeiksin| [2009). It is also crucial for 
understanding the abundance and distribution of heavy elements in 
the IGM, whose ionization state depends sensitively on the local 
metagalactic radiation field (e.g., Song aila 1998. 120051 ; iKim et alj 
l2002bl ; lAguirre et alJl2004 l ;lBolton & ViefcOllh . Additionally, the 
ionizing background is an important input parameter for cosmo- 
logical simulations because it regulates the dominant heating an d 
cooling in the IGM jDave et all 19 99; Spring el & Herng uist 2003), 
which forms the fuel supply for later galaxy formation. Finally, 
the ionizing background holds important clues about galaxies and 
quasars, because they are the dominant sources behind it. Precise 
measurements can constrain the star formation rate, the escape frac- 
tion of ionizing photons from galaxies, and the importance of lu- 
minous quasars (Madau et al. 1999; Faucher-Giguere et al. 2008a, 
120091; lHaardt & Madaull2012h . 

Perhaps most importantly, the ionizing background is tied 
inextricably to the reionization process, when the global ion- 
ization state of intergalactic atoms changes rapidly. For exam- 
ple, measurements of the Hi ionizing background at z ~ 5- 
6 show that hydrogen reionization appears to proceed rela- 
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tively slowly dBolton & Haehnelt 2007). Its properties will also 
be crucial for un derstanding He II reionization, which i s due 
to bright quasars |Sokasian et alj 1200 3l ; iFurlanetto & Oh] 120081; 
IMcOuinn et all [2009). Based on studies of the effective opti- 
cal depth of the Hell Lya forest, the reionizati on of He II in 
the u n iverse seems to h a ve completed at z ~ 3 (Reimers et al. 
1 19971 ; iKriss et alj l200ll ; IZheng etalJ |2004 IShull et alj |2004|) . 
The evolution of the ionizing background during and after 
He II reionization is critical to interpreting new and upcom- 
ing He II Lya fores t results from HST/C OS dShull et alj bold ; 
IWorseck et alj l201ll ; ISvphers et al] l2012h . Theoretical calcula- 
tions have attempted to address this evolution by semi-analytic 



modelling dDixon & Furlanettol 1 20091 ; IFurlanetto & Dixor] 
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and hydrodynamic simulations of the IGM dSokasian et 
iBolton et al.ll200fj ; IPaschos et ai]|2007l ; IMcOuinn et alj|2009 

There is a long history of calculations to e stimate the prop- 
erties of the metagalactic ionizing radiation field. lHaardt & Madaul 
dl99fj|) made a landmark study of the ionizing background using 
a cosmological radiative transfer model for ionizing photons trav- 
eling through a clumpy IGM. By combining state-of-the-art con- 
straints on the distribu tion of ionizing source s and the absorber dis- 
tribution of the IGM, Haardt & Madau dl99fj) were able to compute 
the evolving ionizing background of H I and He II . Further studies 
dFardal et alJI 19981 ; iFaucher-Giguere et alj|2009l ; lHaardt & Madaul 
120121) have updated this framework with new constraints on the 
population of ionizing sources and the distribution and properties 
of IGM absorbers. However, all of these studies treated the ion- 
izing background (and its sources and sinks) as uniform compo- 
nents, which is a reasonable approxim ation for the H I background 
(at least at low and moderate redshifts; Meik sin & White! 20041) but 
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is a poor approximation when bright, rare sources dominate the 
emissivity (as is the case for quasars and the He II ionizing back- 
ground}; 

lFardaletaI] [l998) showed how the relatively large mean sep- 
aration of He II ionizing sources could contribute to the significant 
observed fluctuations in the ionizing background and hence in the 
observable He II Lya effective optical depth. An analytic descrip- 
tio n of variatio ns in the metag alactic radiation fie l d was introduced 
by liuol h992h. expande d by iMeiksin & White! J2003I) . and later 
used by Furlanetto ( 2009) to study fluctuations in the He II ionizing 
background. Despite this theoretical interest, there has been no ef- 
fort to include the effect of these fluctuations on the ionizing con- 
tinuum opacity within a cosmological radiative transfer model. In 
this work, we attempt to show the self-consistent effect of these 
fluctuations on the mean ionizing background. 

We begin in Section 2 with a description of our implemen- 
tation of a cosmological radiative transfer model to calculate self- 
consistently the He II ionization rate. Then, in Section 3, we present 
the results of our model. In Section 4, we use the results from that 
model to calculate the evolution of the He II effective optical depth 
and compare it to observations. We discuss our model assumptions 
and compare to previous work in Section 5. We conclude in Section 
6. 

In our calculations, we assume the following cosmology: 
Cl m = 0.26 Ov = 0.74, Q b = 0.044, and h = 0.74 
dDunklevetaT1 l2009). All distances are given in comoving units 
unless otherwise specified. 



2 INPUTS/METHODS 

2.1 Cosmological Radiative Transfer 

To calculate the He II ionizing backg round, we employ a cosmo- 
logical radiative transfer model lHaardt & Madaulll996h . By con- 
sidering photon conservation in a comoving volume element, the 
specific intensity of ionizing radiation J v behaves as 
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where H(t) is the Hubble parameter, c is the speed of light, a v is 
the absorption coefficient (with dr v — a v dl p and dl p the proper 
line element), and e v is the emissivity. This approach assumes that 
each comoving volume element can be described as an isotropic 
source and sink of radiation through e„ and ct v , respectively: we 
will revisit this assumption later on. The solution to the cosmolog- 
ical radiative transfer equation is 



J„ (zo) = — 



' dl (1 + z f 

dz (1 + z) 3 ex P[ -r ( l 'o, «o, z)\. (2) 



where dl/dz — c/H(z) is the comoving line element, v — Uo(l + 
z)/(l + zq), and f is the effective optical depth experienced by a 
photon at frequency and redshift zq since its emission at redshift 
2. f is calculated using e~ T = (e~ T ) averaging over all lines of 
sight. For Poisson-dis tributed absorber s with Hi column density 
Nu i this opacity is JParesce et al Jl 1980h 
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dz 
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where d 2 N/dNmdz = f(Nni,z) is the column density distri- 
bution function (CDDF) of neutral hydrogen absorbers. The most 



common simple form of the CDDF is a power law in column den- 
sity and redshift: /(Ahi, z) oc (1 + z) 1 , but we will allow 
more sophisticated models as well (see §[2A3j- 

The optical depth of an absorber to ionizing photons of fre- 
quency v is given by 



(4) 



T v = Ah 1(7 Hl(f ) + iVHclOHclM + Ahc II^Hc II (V) , 

where AT, are the column densities and <Ji are the photoionization 
cross-sections of ion i. Because only the column density distribu- 
tion of Nu i has been measured, we use a model for the relationship 
between Ah i and Ahc ii to calculate the He II ionizing opacity (see 
§ 12.1. It . In the frequency range contributing to the Hell ionizing 
background {y > vuaii = 4 x mi) we assume the contribution 
to the optical depth from He I is negligible. Finally, the ionization 
rate for species i is given by 



r i (z) = 4tt 



Jy{z) 

hv 



Oi{v)dv, 



(5) 



where Vi is the ionization threshold of species i. 

Because the ionizing background and the continuum opacity 
are interrelated through the conversion of JVh i to A r H C n described 
in the next section, the procedure must be iterated over the entire 
redshift range until convergence. 



2.1.1 Absorber ionization Structure: Nn i to Nu_ c n 

The relationship between TVh i and A^Heii is usually param eterized 
by the quantity r\ = Nucu/Nhi dMiralda-Escuddl 1993h . In the 
optically thin case, r\ is given by 



Thi o&n Y 

rHell Out iX ' 



(6) 



where i and aj| c n are the case-A recombination coefficients of 
H I and He II respectively. In an optically thin environment, pho- 
tons produced by recombinations to the ground state of He II will 
escape from the local medium, hence our choice of case-A recom- 
bination coefficients. Note, however, that these coefficients enter 
only in the ratio, so this choice does not have any significant effect. 

To more generally translate Hi column densities into Hell , 
we adopt a fit to numerical simulations that accounts for 
self- shielding in neutral hy drogen systems dFardal et alj 1 19981 : 
iFaucher-Giguere et alj|2009n . 



Y 



T"HI 



-/HI = THcII + 



T"HoII 



-ike II. 



(7) 



16X1 + Ami 1 + BTHeii" 

where n = (TiNj, A = 0.15 and B = 0.2 are fitting coefficients 
used by Faucher-Giguere et a D <2009h . and h = Ti/n e OLi with 
n e = 1. 4 x 10 " ^ cm- 3 (A^Hi/10 172 cm- 2 ) 2/3 (r H i/10- 12 
s" 1 ) 2 / 3 JSchavd 1200 ll) . At small Hi column densities 
(iV H i<10 15 cm -2 ), A'hc ii = ^thm^Hi as expected. 
He II becomes optically thick to ionizing radiation for larger 
column densities (A^hi ~ 10 15 -10 17 cm -2 ), so r\ increases by 
a factor of a few as more Hell forms while hydrogen remains 
highly ionized. Then, for Ahi^IO 17 cm -2 , r\ steeply drops as 
the systems become optically thick to H I ionizing photons. 

For systems with AThi > 10 ls cm -2 , the numerical fit sys- 
tematically under-pred icts the amount of He II from the original 
model (see Figure 1 of IFaucher-Giguere et alj|2009h . For frequen- 
cies near fHeii, the opacity is unaffected because these high TVhi 
systems are still optically thick thanks to H I absorption. However, 
for i/>2.5^Hcii, absorbers with AThi 
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start 
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to become optically thin due to their relative lack of He II . Fortu- 
nately, the total ionization rate only changes slightly because the 
range of affected column densities is small and the vast majority of 
ioni zations occur at lower fre quencies (cthcII oc v~ 3 ). 

lHaardt & Madaul d2012h applied a similar method to fit the 
absorber structure that considers the average V within absorbers 
instead of the external "optically-thin" V. While their method pro- 
vides a better fit to the numerical models at A^ hi > 10 18 cm -2 , 
it differs from the Faucher-Giguere et ai] j2009t) model only in the 
details for the more important thoII ~ 1 (iVki ~ 10 16 cm -2 ) ab- 
sorbers. This is an example of one of the systematic uncertainties 
in our procedure: these models for rj must assume physical char- 
acteristics for the absorbers (densities, temperatures, and geome- 
try, for example) that are both uncertain and simplifications of the 
true IGM physics. For concretenes s, the numerical absorber model 
from lFaucher-Giguere et al" (2009) assumes uniform density semi- 
infinite slabs with a thickness determined by the local Jeans length 
(at T = 20, 000 K) in photoionization equilibrium with both an ex- 
ternal radiation background and internal recombination processes. 



2.1.2 Recombination Emissivity 

Recombinations of He III to the ground state of He II will produce 
ionizing continuum radiation. Although the recombination rate in 
a uniform density medium can easily be estimated from ionization 
equilibrium, the real universe requires a more detailed approach 
for two reasons. First, density inhomogeneities in the IGM sub- 
stantially boost the recombination rate. We can model this by inte- 
grating over the H I column density distribution of the Lya forest. 
Second, recombination photons produced inside optically thick ab- 
sorbers will not escape to affect the IGM. 

We model the recombination emissivity of IGM ab- 
so rbers with a numeric a l fit to the radiative transfer models 
of iFaucher-Giguere et ID (§009). The emergent specific intensity 
from an absorber with He II column density 7Vh c ii can be approx- 
imated by 



/r^Hell) 
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l HeII, 



(8) 



where the second factor is the fraction of ionizations to the ground 
state and the local ionization rate is Thou- Nt = 10 17 ' 3 cm -2 
is the approximate threshold He II column density above which the 
emission becomes saturated by absorption within the absorber itself 
(the decline at larger columns is approximated by the last factor). 
4>v,icc is the normalized recombination emission profile: 

4>u,iec cx v~ x e~ h "' hBT 0{v - men), (9) 

where 6(x) is the Heaviside step function. The effective frequency 
width of this emission is Av/v ~ ksT //ifHell ~ 0.03, which 
limits the distance these photons can travel to < 30 Mpc before 
redshifting below the He II ionizing edge. 

The total comoving emissivity from recombinations is then 

e^recfz) = -jr(l + z)~ 2 
at 

x / dNnif(Nni,z)ATrI I i ; c (Nnen), (10) 
Jo 

where the intensity depends implicitly on Nh i through the conver- 
sion factor r\. 

We note here that the recombination photons can have a much 



larger effect on the ionizing background than one might naively 
expect from their emissivity. As we shall see later, increasing the 
emissivity also increases the mean free path of ionizing photons, 
which amplifies the effect of the additional ionizing photons. We 
will explore this issue further in § 13.1.51 



2.2 Mean Free Path 

The opacity per unit redshift, df/dz, was integrated in equation l[3} 
to calculate the total opacity between two redshifts: 



df 

dz 



dN HI f(N H i,z)(l-e- T "). 



(11) 



At a given redshift, df/dz describes the local opacity due to the 
forest of individual absorbers in the IGM. By inverting this quantity 
and converting from redshift difference to a comoving distance, we 
find the distance per unit optical depth, which is simply the mean 
free path: 



W",*) = £(f) ■ 

If f(Ni,z) = N N7 B (l + z)-> and a* 
ing mean free path reduces to 



cro{v/vi) 
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, the comov- 



A m f p (f, z) 



09-1) 



r G (2-/3)jV, 

i 



3(0-1) 



(13) 



where To is the Gamma function. The redshift dependence of the 
mean free path in this simplified model is then A m f p oc (14- 
z) - (t+ 1 - 5 ). This power law dependence is a good approximation 
to describe the evolution of the mean free path of H I ionizing pho- 
tons (Ah i) in our model because the H I CDDF is fixed, but we find 
that it fails to capture the more complicated Thc ii-dependent evo- 
lution of the mean free path of Hell ionizing photons (AhcIi; see 

Recent efforts by IProchaska et alj j2009h and lO'Meara et alj 
d2012l) have directly measured the Hi ionizing mean free path of 
the IGM near z ~ 4 and z ~ 2 respectively. For an identical 
distribution of absorbers, they would report different values than 
obtained by our approach because they define the mean free path 
as the distance traveled by a photon through the evolving IGM and 
redshifting as it does so, rather than the path that could be traveled 
if the IGM and photon retained their original properties (as is the 
usual definition for theoretical work). We follow the latter defini- 
tion here. 



2.3 Fluctuations 

In a smooth, fully-ionized IGM, the intensity of ionizing radia- 
tion from an individual quasar falls as exp[— r/A m f p ]/r 2 . Given 
a distribution of quasar luminosities and a mean free path, a prob- 
ability distribution of intensities can be computed ass uming ran- 
dom placement of quasa rs following Poisson statistics feud 19921 : 
Meiksin & White 2003). The effects of this distribution on the 
mean ionizing background have not previously been considered. 
The next stage in our model is therefore to incorporate the distribu- 
tion (in a somewhat ad hoc manner) in order to understand better 
the implications of this fluctuating background. 

We use the Hopkins et al. (2007) B-band quasar luminosity 
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Figure 1. Distribution of ionization rates at z = 3.0 for A m f p 
30, 60, 120 Mpc (long-dashed, dashed, and solid, respectively). 



function to describe the distribution of relative quasar luminosi- 
ties, assuming an average quasar spectral energy distribution such 
that the specific luminosity at the H I ionizing edge is proportional 
to the B-band specific luminosity {Lb), then extrapolating to the 
He II ionizing edge by a spectral index a. Additionally, we con- 
volve the quasar luminosity function with a distribution of far- 
ultraviolet spectra l indices that roughly matches observations by 
iTelferetalJ 120021) : a Gamma distribution over 0.5 < a < 3.5 
with a — 1.5 and a a — 0.5. 

We use th e method of characteristic functions from 
iMeiksin & White] d2003h to determine the probability distribution 
of intensity, f(J ), then scale linearly to T by F = J x (r)/(J) 
dFurlanettofe oOO). The last assumption of proportionality between 
the intensity of radiation and the ionization rate is not strictly true; 
the intensity at higher frequencies should be more uniform because 
the mean f ree path is much larger, al though the effect is modest 
in practice jDixon & Furlanettd F2012). In our calculation of /(F) 
we use the mean free path of the "average" He II ionizing photon, 
AhgII = A m f p (z/), where v is defined by 



v x — <rnen(y)dv, 

hv 



(14) 



in an attempt to average over the frequency dependence of the back- 
ground fluctuations. In general, AHeii is substantially larger than 
AhgII, so our approach provides a conservative estimate when used 
to calculate the amplitude of ionizing background fluctuations. 

Figure Q] shows how the f(T) distribution varies with mean 
free path. When the mean free path decreases, the peak of the dis- 
tribution skews towards smaller F relative to the mean. For T below 
the mean, the Hell opacity of each absorber will increase. Because 
most of the Universe is in this configuration - while only a small 
fraction experiences a very large ionizing background - we expect 
the mean T to decrease. We explore this effect in the following dis- 
cussion. 

We incorporate these fluctuations into our ionizing back- 



ground model by averaging the opacity dr / dz (equations [3] fTTT l 
over the distribution f(T): 



dz I 



■ lf \ = 1 £/(r)dr, 



dz ' 



(15) 



where df/dz implicitly depends on T and f(T) is initialized with 
the mean free path calculated in the uniform model. This process is 
repeated using the same f(T) for each frequency in equations (O 
and (O and separately at each redshift to calculate a new Fhb ii and 
iterated until convergence. 

Unfortunately, as presented above, the Fh c ii calculation does 
not converge. At relatively high redshifts (z > 3.5) the mean free 
path is short enough (A m f p < 20 Mpc) that integrating over f(T) 
greatly increases the opacity. This increased opacity at high red- 
shift propagates small values of Tn e 11 to lower redshifts, and the 
iterative effect pulls T down to zero at all redshifts. 

The reason for this apparent divergence is actually obvious: 
our cosmological radiative transfer model assumes that ionizing 
photons are emitted uniformly throughout the universe (with a con- 
stant e„ in equation QJ, but the real quasar sources are of course 
point-like. Since the ionizing background near a source is much 
stronger than the average, the local IGM will be less opaque to 
ionizing photons. This increased transparency near the emitting 
sources allows ionizing photons to travel an average distance, 
A m in > AhcII, before being absorbed, which should not be subject 
to averaging over regions that are devoid of sources. In this discrete 
source limit, we should approach the distribution of the ionizing 
background differently as well. Qualitatively, we can account for 
this picture by considering only these regions around the sources 
when propagating ionizing photons - the more opaque regions ex- 
terior to these "bubbles" have very few ionizing photons, so their 
properties are not relevant to the average. 

To implement this picture, we set A m i n (see below for details 
on its calculation) to be the smallest allowed mean free path when 
generating f(T). Additionally, we modify equation dl5t so as to 
consider only those T consistent with AhcII > A m i n : 



df \ 
dz) 



df 
dz 



f(T)dT + 



(16) 



where r m i n is defined by AHeii(r m i n ) = A m i n . Larger values of 
F correspond to regions inside the illuminated bubbles so are still 
allowed. 

Once the opacity due to the fluctuations in the background 
dominates, T m i n ~ f , and the first term in equation l l 1 6t will domi- 
nate the average opacity. In this regime, the resulting ionizing back- 
ground is constructed from A m i n -sized bubbles around each quasar. 
However, (by ignoring regions outside the zones near quasars) the 
method implicitly assumes that these near-zones fill space, which is 
not the case. The filling factor of these A m i n -sized bubbles is given 
by 



fx 



(17) 



where Z/ m in is our assumed minimum quasar luminosity; we take 

-L'min — 

10 43 erg s" 1 for o ur calculati ons, but the results are not 
sensitive to that choice; see iFurlanettd j2009h . fx min declines be- 
low unity for z > 2 (~ 0.4 at z ~ 3), implying that we overesti- 
mate the true background. On the other hand, we can obtain a lower 
limit for the background by assuming that Thc ii is zero outside of 
these regions, which amounts to multiplying our modelled value 
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Figure 2. Left: C olumn density distribution fun ctions f{N H 1,2 = 2.5) considere d in the text: iRudie et all l2012al) (solid black ), Haardt & Madaul J2012b 
(dot-dashed blue), Faucher-Giguere et al. (2009) (dotted red), O'Meara et al 1 120121) (long-dashed green). IWorseck & Prochaskal 1201 ll) (short-dashed-long- 
dashed purple). The vertical dashed line shows the Nm corresponding to a Hell "LLS". Right: Relative contribution to the continuum opacity at ^Hcll P er 
log(AT H i). 



by the filling factor /A min ■ We will show both alternatives as part 
of our main results. 

To determine A m i n , the average distance an ionizing pho- 
ton (emitted from a q uasar) can travel before being absorbed, 
we adopt the method of iFurlanetto & Oh (2008). We assume that 
the gas density probability distribution -P(A) is described by the 
Miral da-Escude et al.1 d2000h model and that all absorber systems 
with d ensity greater than a threshold density Ai remain Hell . Fol- 
lowing [MiraldajEscudiletj^ ( feOOOl) . we define the mean free path 
of ionizing photons as the average distance between regions with 
A > A;. The distance between these regions is then 



A, 



Ao[l — Fy(Ai 



-2/3 



(18) 



where Ao is a normalization factor, and F v (Ai) is the volume fill- 
ing factor of ionized (A < Ai) gas. The power law index, —2/3, 
comes from assuming spherical absorbers. The average distance 
an ionizing photon travels from its source quasar, A m i n , is then 
the point where Ai(Ai) equals the distance from the quasar at 
which all densities below Ai are ionized. We find for our assumed 
quasar SED and quasar luminosity function a luminosity-weighted 
(A m in) ~ 35 Mpc, which we apply to the fluctuating background 
prescription described above. 



2.4 Model Input Parameters 

Other than our simple model assumptions, the largest sources 
of uncertainty in our analysis are three observed parameters: the 
He II ionizing emissivity (e„), the H I ionization rate (Th i), and the 
neutral hydrogen column density distribution (/(TVhi, z)). In this 
section, we discuss the range of observed values for these parame- 
ters. 



2.4.1 Hell Ionizing Emissivity 

We adopt the Lyman -limit quasar ionizing emissivity from 
lHaardt & M adau (20ll|), 



£912(2) = 10 ' ergs Mpc Hz 

.4.68 exp[-0.28z] 



exp[l. 77z] +26.3' 



(19) 



w hich is a fit to the int egrated B-band quasar luminosity function 
of iHopkins et al.l ( 2007) converted to 1 by a constant factor, 



L.hi =Lb x lO^ergs^Hz" 1 ( ^ 

\ Lb 



(20) 



This factor is effectively an estimate of the average quasar spectrum 
between vb and i/hi. For frequencies above the Lyman limit, we 
assume a power law spectrum with e v oc v~ a . For reference, the 
integrated quasar emissivity given by equation ( |19t increases by 
- 30% from z = 3-2. 

The uncertainty in the He II ionizing emissiv ity is a combina- 
tion of the uncertainty in the lHopkins et al. U2007I) quasar luminos- 
ity function and the assumed average quasar spectrum. The former 
uncertainty is likely to be small, because the integrated quasar B- 
band emissivity at z > 2 comes pred ominantly from the brightest, 
and therefore best measured, sources dHopkins et alj20071) . The lat- 
ter unc ertainty is domina ted by the choice of far-UV spectral in- 
dex a. lTelfer etaT] ( l2002l) find a = 1.57 ± 0.17 for a composite 
spectrum of 77 radio-quiet quasars, while the composite includ- 
ing an addi tional 107 radio -loud quasars has a — 1.76 ± 0.12. 
In contrast, IScott et al.l (2004J) found that the average spectral in- 
dex for the ir sample of 85 so urces was considerably harder, a = 
0.56 + q jg. lShull et alj J2012t> measured a best-fit spectral index of 
q = 1.41 ± 0.21 for their sample of 22 sources using HST/COS. 

We adopt a — 1.6 as our fiducial value. Note that, because the 
Hell Lyman limit ^Hcii = 4i^hi, a change in the spectral index 
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Aa corresponds to a factor of 4 

^HcII- 



2.4.2 Hi ionization Rate 



difference in the emissivity at 



The absorber model in § 12. 1.1 1 depends on the Hi ionization rate, 
Thi- Measurements of Thi from z ~ 2-3 yield val ues ~ 0.5 — 
1.0 x 10~ 12 s -1 from flu x decrement observations dRauch et alJ 
19971: iBolton et al.1 I200I; [McDonald & Miralda-Escudel l200ll ; 

i> o 



Faucher-Giguere et alj|2008j > or ~ 1.0 - 3.0 x 10 _i/ s" 1 from 
proximity effect measurements {Scott et al. 2000) . The mos t recen t 
cosmological radiative transfer model by Haardt & Madaul d2012h 
suggests Thi ~ 0.8 — 0.9 x 10" 12 s _1 , but as discussed in the 
next section, they may have significantly underestimated the total 
Hi opacity of the IGM. W e adopt Thi = 0.5 x 10~ 12 s" 1 , the 
value measured recently by Faucher-Gigue re et al. (2008a), as our 
fiducial value but consider a range of plausible values. 



2.4.3 Column Density Distribution 

The column density distribution of neutral hydrogen /(Ah 1, z) = 
d 2 N/dNuidz has been measured several times and over a range 
of redshifts through observations of the H I Lya forest. Early ob- 
servations indicated that the A?hi distribution is well-fit by a 
power law of the form f(Nui) oc A^f with f3 ~ 1.5 over a 
wide range of observed column densi ties (10 12 < A?hi < 10 22 
cm" 2 ) and redshifts (z ~ 0.2-3.5) jTvtleil 1 19871) . Recent stud- 
ies of H I ionizing continuum opacity in stacked quasar spectra 
at z ~ 2 and z ~ 4 suggest a deficit of Lyman-limit sys- 
tems (10 17 ' 2 < Nhi < 10 19 cm -2 ; LLS) and intermediate 
Hi column density systems (10 15 < JVhi < 10 17 ' 2 cm -2 ) rel- 
ative to the canonical single power law model, and several authors 
have pr oposed multi-step power law distributions to describe this 
feature dProchaska eTal. 2009, 2010; Worse ck & Prochaska|[2oTll ; 
lO'Meara et afll2012l ; lHaardt & Madaill2012h . iRudie et alj d2012ah 
performed the largest survey of 10 12 < iVm < 10 17 ' 2 cm -2 
systems to date and found no evidence of the deficit suggested 
by stacked quasar spectra studies. They found that their measured 
distribution is well-parameterized by a relatively steep /3 ~ 1.66 



power law for Nni < 10 cm" and a ~ 1.48 power law for 
larger Hi column densities. The left panel of Figure [2] shows sev- 
eral of these distributions graphically. 

The redshift evolution of the CDDF is usually parameterized 
by a power law f(Nui,z) oc (1 + z) 7 . However, observation- 
ally this 7 appears to depend on ATht, implying that the shape 
of the CDDF evolves with time. The observational constraints on 
7 for z > 2 in the Lya forest reg ime (A^hi < 10 1 7 2 cm -2 ) are 
7 ~ 2.0-3.0 from line-counting dKim et al 2002ah and measure- 
ments of the effective optical depth ( Faucher-Gig uere et all 2008b; 
iDall' Aglio et alj2008l) . The number densities of super-Lyman limit 
(10 19 cm " 2 < ATri < 10 20 - 3 cm _2 ) and damped Lya (JV HJ > 
10 2 3 cm ~ 2 ) absorbers appear to evolve more slow ly with 7 ~ 
1.7 dO'Meara et all200llWorseck & Prochaskj20i"il) and ~ 1.27 
dRao et al.l2 006). respectively. IRudie et alj l2012ah found that their 
data were consistent with 7 = 2.5 and 1.0 for A^hi below and 
abo ve ~ 10 15 cm" 2 , respec t ively. 

IWorseck & Prochask j d201ll) and lHaardt & Madlul d2012h 
compiled these observations and constructed similar multi-step 
power law CDDFs. The primary difference between the two is 
the enhanced redshif t evolution (t_ = 3.0) of Lya forest ab- 
sorbers in the lHaardt & Madaul d2012h model compared to the 
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Figure 3. Uniform and fluctuating Thc ii in the fiducial model (solid and 
dashed curves, respectively), the minimum average ionization rate in the 
"plateau" regime rjj c ii x f\ (dot-dashed), and the contribution from 
recombination emission (long-dashed). The shaded region represents values 
of Tjie 11 that are "reasonable" in the context of our near-zone model. 



IWorseck & Prochaskal d201lh model (7 = 2.04). Both models de- 
termine the redshift evolution of the CDDF by comparing to ob- 
servations of the evolution of the Hi Lya effect ive optical depth, 
which is proportional to (1 + z) 7+1 . However, Haardt & Mada 
20121 ) calibrat e to the measurement s of [ Faucher-Giguere et 



dau 
"aH 



2008b]), while IWorseck & Prochaskl J20T J) chose the measure- 



ments of IDall' Aglio et al ] d2008h . It is unclear why such a differ- 
ence exists in the effective optical depth evolution measured by 
these two groups, but it does not significantly affect our results. 

In the following sec tions, we use the broken power-law CDDF 
from lRudie et al. (2012a) as our fiducial model. Their model repre- 
sents the first solid measurement of intermediate H I column den- 
sity absorbers that are critical to the He II ionizing opacity, and 
is consistent with measurements of the Hi Lya effective optical 
depth (G. Rudie, priv. comm.). However, as the following sections 
will show, our choice of CDDF does not have significant implica- 
tions for our final results, given the overall uncertainty in the am- 
plitude of the ionizing background. 



3 EVOLUTION OF THE IONIZING BACKGROUND 
3.1 The Ionizing Background With Uniform Emissivity 

The solid curve in Figure [3] shows how the He II ionization rate 
(Tire 11) evolves in our fiducial model, ignoring fluctuations in the 
ionizing background. The uniform background model results in a 
steeply evolving ionizing background from z ~ 3-2, with an ion- 
ization rate that increases by a factor of ~ 2 over that range before 
flattening out substantially at later times. In the following sections, 
we discuss how variations in the input parameters affect this result. 
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Figure 4. Top: Tfjeii in the uniform background model calculated for 
CDDFs from Figure [2] Bottom: Effect of assumed average quasar spec- 
trum shortward of 912 A , given by F v oc v~ a (solid curves), and assumed 
(constant) Th i (dot-dashed curves). 



3.1.1 Column Density Distribution 

We considered a variety of CDDFs in our model. The top panel 
of Figure [4] shows the uniform ionizing background calculated 
with CDDFs used in ea rlier ionizing background calcula tions by 
lHaardt & Madaul d2012l) and lFauc her-Gigu ere et"aj] d2009h . the di- 
rect measurement at (z) ~ 2.4 by iRudie et alj d2012al) . and in- 
direct extrapolation from higher redshift opacit y measurements 
dWorseck & Prochaskj|201 ll;lo'Meara et alj|2012h . In general, de- 
spite the significant differences between CDDFs apparent in the left 
panel of Figure [2] the evolution of the uniform background from 
z = 3-2 is fairly insensitive to the CDDF The most significant dif- 
ferences are due to the different redshift evolution of the CDDFs, 
which is not very well constrained. 

The right panel of Figure [2] shows the relative contribution 
to the ionizing continuum opacity at the Hell edge as a func- 
tion of Nhi- Most of the opacity is due to Hell "LLSs" with 
A^Heii ~ a HeW rjut me multi-step power law CDDFs have an in- 
creased contribution from relatively low A^ hi ( 5 10 15 cm" 2 ) ab- 
sorber s compared to the shallow power law Faucher-Giguere et al. 
(2009) CDDF. The Hi column density corresponding to the peak 
Hell opacity contribution varies from 10 15 to 10 16 - 7 cm de- 
pending on the shape of the CDDF. 

Figure [4] also shows that the normalization of ThcII depends 
sensitively on the total opacity calculated from the CDDF, which 
can vary significantly between models. If THeii were accurately 
measured near z ~ 2, that measurement could in principle be used 
to help distinguish between models. However, measuring THeii 
directly is extremely difficult, and as shown in the following sec- 
tions, the other model parameters can be adjusted to produce simi- 
lar differences in the normalization. For example, measurements of 
both the 7/thin parameter and Th i could potentially be used to con- 
strain acceptable normalizations of FhcII (because the expected 
He II Lya opacity in the IGM depends strongly on the value of the 
former parameter; see equation |6), but the current constraints on 



these parameters are too weak, and the degeneracies are too strong, 
to distinguish between the models presented in this and the follow- 
ing sections. 

3.1.2 Quasar Spectrum 

The solid curves in the bottom panel of Figure [4] show how the 
range of observed values of the quasar spectral index a affects 
the Hell ionization rate. A harder spectrum, which produces more 
ionizing photons at ^Hcii, results in a higher ionization rate. Fix- 
ing the emissivity at z^Hcii and changing the spectral index has 
very little effect on the resulting ThcH- In contrast, we find that 
ThcII changes more strongly than linearly with eHcii; this is 
because the absorber structure changes with the ionizing back- 
ground (and hence the emissivity). In general, as THe n increases, 
the Hi column density corresponding to a Hell LLS increases. 
Since Ahi/(Ahi, z) is a decreasing function of Ahi, the num- 
ber density of He II LLSs, and thus the overall opacity, decreases. 
This behaviour is sim ilar to the emissivity-r feedback studied by 
iMcOuinn et alj d201 ll) . The redshift evolution of the background is 
affected by a as well, but the effect is subtle. 

3.1.3 H I ionization Rate 

The dot-dashed curves in the bottom panel of Figure [4] show how 
the Hell ionization rate is affected by the assumed value of Fhi. 
The effect is similar to changing the number of He II ionizing pho- 
tons, because both parameters modulate the ratio of He II to H I in 
absorbers. While the decrease in Hell opacity with an increas- 
ing number of He II ionizing photons is straightforward in princi- 
ple, the relationship between Th i and Thc ii is more subtle. Con- 
sider an optically thin absorber: if Th i decreases, the amount of 
Hi in a fixed physical structure will increase while the amount 
of Hell stays the same. This shift of the Hi column density cor- 
responding to a He II LLS causes Fh c ii to change inversely with 
Thi; that is, if Thi is larger, the calculated ThoII will be smaller, 
and vice- versa. Thi appears to affect the redshift evolution more 
strongly than a. 

3.1.4 Mean Free Path 

The solid curve in Figure [5] shows the evolution of AHeii in the 
uniform model. We also show how Ah i increases with cosmic time 
(dotted curve); for ease of comparison we scale this curve to AhcII 
at z — 2. In contrast to the power-law evolution of Ah i (described 
by equation \\3[, XhbII evolves much faster than a simple power 
law. 

The evolution of the mean free path at the He II ionizing edge 
in our fiducial model is well-approximated by a power law with an 
index that itself evolves as a power law, 

AHeii ~ 210 comoving Mpc x — ■ — (21) 

C(z) = -2.29 x (-^-J ■ (22) 

This fit differs by no more than ~ 3% from our full numerical cal- 
culations over the redshift range z = 2-3.8, but we caution the 
reader that the systematic variations from our model input param- 
eters are much, much larger than this. We also caution the reader 
against using this fit at z > 3.4, where fluctuations in the ionizing 
background must be included. 
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due to an inhomogeneous IGM through the so-called clumping fac- 
tor, C = (rune) I ({rii) {n e )) . Usually, this is estimated from sim- 
ple phenomenological arguments or from the density structure in 
numerical simulations. However, these approaches are not entirely 
satisfactory, as the clumping factor should incorporate information 
that depends on the distribution of ionized and neutral patches. For 
example, recombinations that occur inside of dense, self-shielded 
systems do not produce photons that can ionize the IGM, as the 
resulting photons are trapped within the systems. 

With our detailed model, we can estimate this factor for 
He III self-consistently (given a model for the emitting and absorb- 
ing populations) by explicitly following the fraction of recombina- 
tions that occur inside of self-shielded systems. In particular, we 
have 



z / (hv)dv 



I )(nHeIIl)(»1e) ' 



(23) 



which describes the effective recombination rate after correcting 
for self-absorption of ionizing recombination photons within the 
emitting clouds relative to a uniform IGM. In our fiducial model, 
C c a increases from C e g ~ 1.0 at z = 3.4 to C e g ~ 3.7 at z = 2. 



Figure 5. Evolution of the Hell ionizing photon mean free path with 
redshift (black), evaluated at the average ionizing energy, for the uniform 
(solid) and fluctuating (dashed) models. The hydrogen ionizing photon 
mean free path is shown as the dotted curve, scaled to the He II mean free 
path at z = 2. The red long-dashed curve shows the average separation be- 
tween luminous (vbLb > 10 Lq) quasars given by the lHopkins et alj 
120071) OLF. 

Because AhgII is linked to THeH through the absorber struc- 
ture prescription, it evolves more quickly than Ah I- That is, increas- 
ing the mean free path increases the ionizing background, which 
will then increase the H I column density at which He II becomes 
optically thick, which in turn increases the mean free path, etc. This 
feedback effect is the fundamental source of the rapid evolution we 
see in THeH- (In fact, one could argue tha t it is strange that su ch 
rapid evolution does not occur in Th i; see lMcOuinn et alj201ll .) 

The dependence of the mean free path on frequency is a 
function of the logarithmic slope of the CDDF, A m f p oc u :if/3 ~^ 
(equation 1 1 3b . The Hell CDDF is not precisely defined in our 
model, but a mapping of our fiducial Hi CDDF through our ab- 
sorber prescription results in /3hcII ~ 1-43 for the absorbers that 
contribute the bulk of the opacity near the He II ionizing edge 
(10 14 ' 5 < Nhi < 10 17 - ctrT 2 as in Figure [2}, and consequently 
AhcII oc u 1 ' 3 for 1 < v/vneii ^2.5. P/i/Heii ~ 1-35 is typical 
for our fiducial model, so AHell/AHell ~ 1.48. 



3.1.5 Recombination Photons 

Figure[3]shows that the fractional contribution of recombinations to 
ThcH is fairly minor, decreasing from ~ 0.12 at z = 3 to ~ 0.07 
at z = 2. However, because the absorber population is sensitive to 
the emissivity, the fractional increase in F ne n from including re- 
combi nations is larger (~ 0.4-0.2). While iFaucher-Giguere et all 
(2009) found that including recombinations increased Thgii by 
only ~ 10%, Figure [2] shows that their CDDF has a significant 
deficit of the optically thin (Nhi ^ 10 16 cm" 2 ) systems that con- 
tribute most of the recombination emissivity. 

In simple models of the reionization process, it is conventional 
to describe the enhanced recombination rate of ionized species m 



3.2 The Ionizing Background Including Fluctuations 

It is instructive to compare the mean free path from the preced- 
ing section to the average separation between the primary sources 
of ionizing photons, bright quasars with vbLb > W 11 Lq. We 
cal culate the number de nsity of the bright quasars by integrating 
the lHopkins et al.1 d2007t) luminosity function over this luminosity 
range and estimating their average separation by (R) ~ n~ 1//3 . 
The long-dashed red curve in Figure [5] shows this separation; 
(R) ~ 45 Mpc is a good approximation for the entire redshift in- 
terval from z ~ 2-3. 

When the mean free path is similar to the average source sepa- 
ration, fluctuations in the background contribute a substantial opac- 
ity excess. The dashed curve in Figure [5] shows the effect of these 
fluctuations on the ionizing background. Figure[6]shows that, com- 
pared to the uniform model, the fluctuating background model ex- 
hibits a ~ 30-40% dip near Zf ~ 3.1-3.3 for our fiducial input 
parameters and various CDDFs from § 13.1.11 The ionizing back- 
ground for CDDFs with steeper slopes in the He II LLS regime is 
more sensitive to fluctuations and therefore experiences a deeper 
dip. The evolution of all the C DDF models, with the except ion 
of the shallow slope model from lFaucher-Giguere et all I J2009h . is 
very similar. The "feedback" effect between the opacity and the 
ionizing background is weaker for shallower CDDF slopes (e.g. 
iMcOuinn et1ifll201 if), so the net effect of fluctuations is smaller. 

We note that there are two related sources for the differences 
between the curves in Figure[6] the shape of the column density dis- 
tributions and the normalization of the ionizing background. Leav- 
ing the shape fixed, different choices for the B-band emissivity and 
spectral slope result in very similar shapes to those in Figure [6] 
though the location of Zf shifts to lower redshifts as the effective 
background increases. In other words, a higher normalization for 
the ionizing background increases the mean free path and so makes 
the fluctuations less important (except at very high redshifts). The 
more subtle differences in the shapes of the curves are due to vari- 
ations in the shape of the CDDF; this is most dramatically seen by 
the dotted curve (which only reaches the plateau at z > 4.0). 

The most important effect of including fluctuations is to in- 
duce a much more rapid increase in Th c ii with cosmic time. Con- 
sider a region with a smaller than average emissivity. In that re- 
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Figure 6. Top: Fluctuating model ionization rate for the same CDDFs as 
Figure [2] Bottom: Ratio of the fluctuating and uniform model ionization 
rates. 



gion, the ionizing background will also be smaller, so each absorber 
will be more optically thick and the mean free path will be smaller. 
This will further decrease the ionizing background, etc. In a real- 
istic model of the distribution of Thc ii in the presence of quasars, 
most of the volume of the universe has an ionizing background a 
few times below the universal average (to compensate for the very 
brightly illuminated, but small, regions around quasars; see Fig- 
ure!]}, so the mean ionizing background declines. 

The turndown from the uniform model is thus a straightfor- 
ward and robust prediction of our model, though its magnitude de- 
pends on the CDDF. At still higher redshifts, our model shows a 
"plateau" at z > Zf. This flat evolution occurs because the ioniz- 
ing background is dominated by regions around individual quasars, 
which we have described through our A m i n prescription. Because 
the ionizing background within A m i n from a quasar does not evolve 
rapidly over the range z ~ 2-4, in this regime our prescription for 
Thc 11 also does not vary significantly. 

However, as we have emphasized in i]2.3l our default prescrip- 
tion assumes that the mean free path is at least as large as that near 
quasars everywhere in the universe, while in fact the regions around 
quasars fill only a fraction of space 0.5 at 2 > 2.7). The dot- 
ted curve in Figure [3]- our default approach - therefore provides 
a conservative upper limit to the ionizing background in this early 
regime, while the dot-dashed curve - which sets the ionizing back- 
ground in regions outside of these near-zones to zero - sets a lower 
limit. 

These calculations show that the ionizing background can 
evolve very rapidly near z ~ 3, even without any assumptions 
about an evolving Hell fraction. The precise degree of evolu- 
tion is uncertain, but it is at least a factor of a few (even in the 
uniform emissivity model) and likely nearly an order of magni- 
tude when fluctuations are included. In other words, even without 
He II reionization, we should see a rapid increase in the intensity 
of the metagalactic radiation field. We will consider the observable 
implications of this conclusion in the following section. 

An important caveat is that our estimates for ThcH in this 



"plateau" region depend upon the value of A m i n . Our A m i n assump- 
tion effectively sets a "floor" in the ionizing background defined by 
the near-zone properties of individual quasars. This allows ionizing 
photons to travel to reasonable distances from their sources (and 
hence to later cosmic times or smaller redshifts) more easily, which 
is required to sustain a non-zero ThcH when the opacity is subject 
to the contribution of fluctuations. 

However, we caution the reader that such a "plateau" in the 
ionization rate did not likely occur in reality; at early enough times, 
while He II reionization is ongoing, regions outside of quasar near- 
zones will be out of ionization equilibrium and may have a non- 
negligible Hell fraction. Our assumption that the opacity is con- 
structed from individual absorbers will fail in the regime where 
the mean IGM contributes substantial opacity from a significant 
He II fraction, and calculations of the evolving He II fraction are 
required to track the true ionizing background. 

In the post-plateau regime (z < 3.1 in our fiducial model), the 
mean free path at the He II ionizing edge is well-characterized by a 
similar power law within a power law as the uniform model (equa- 
tion[2Tb. 



Arieii ~ 201 comovingMpc x 

/1 + Z\ 2 ' 43 

C(*) = -2.48 x (i±£) . 



l + z\^ 



(24) 



(25) 



The primary difference between the uniform and fluctuating back- 
ground fits is the larger power law index of C( z )> a consequence of 
faster ionizing background evolution. 



4 EFFECTIVE OPTICAL DEPTH 

To gauge the observable import of our results, we will briefly con- 
sider how they manifest in the evolution of the IGM opacity to 
far-ultraviolet photons. Hell Lya ab sorption has been measure d in 
far-ultraviolet spec tra from z ~ 2-4 jDixon & Furlanettol20 09 and 
references therein; IWorseck et al 1201 ll; ISyphers et al .l201lLl2012l) . 
We will compare to the most basic observable from the resulting 
forest of observed absorption features, the average optical depth 
T e ff for the He II Lya transition. We use two different methods to 
estimate T e g; a semi-analy tic model using a gas dens ity probabil- 
ity distribution P( A) as in lDixon &Furlanetto(2009), and a direct 
integration of the He II Lya opacity from the H I CDDF and our 
absorber structure prescription. 

Under the assumptions of a highly-ionized universe in ion- 
ization equilibrium, line opacity dominated by zero-width opti- 
cally thin absorbers, and a power-law temperature-density relation 
T = To A u , the He II Gunn-Pete rson optical depth can be ex- 
pressed as jDixon & Furlanettol2009h 



T"GP 



13.6 



0.142 



To 
10 4 K 

l+2\»/ 2 



n b h 2 

0.0241 



.2-0.7(7^-1) 



(26) 



This "fluctua t ing Gunn-Peterson approximation" (FGPA; 
Weinberg et al. 1997) relates the continuum optical depth to 
the local overdensity A and ionizing background FHeii- The 
systematic uncertainty in tqp due to the above simplifications 
is collapsed into a normalization constant k, which we adjust 
when comparing to observations. We assume an isothermal 
temperature-density relation (7^ = 1) for simplicity, but this 
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Figure 7. Effective optical depth for the uniform (solid), fluctuating 
(dashed), "plateau" (dotted), and "minimum" (dash-dotted) models, using 
the FGPA (black) and CDDF (purple) methods. The black squares are the 
compiled T en - value s fromlDixon & Furl anettol i2009i) , the red triangles are 
re-binned data from Worseck et aljl201 ll) with estimated error bars, and the 
blue pentagons are from lShull et alj 



does not affect our results signifi cantly. Assuming the gas d ensity 
probability distribution given bv lMiralda-Escude et al . (2000), wc 
then calculate r e g by integrating over the density and ionization 
rate distributions: 



dr/(r) 



dAe 



-T GP (r,A) 



P(A). 



(27) 



We normalized the results for the uniform and fluctuating models 
to produce an optical depth of r = 0.5 at z — 2 to roughly match 
observations. These normalizations require k — 1.12 and k = 1.02 
(equation [26} for the uniform and fluctuating models, respectively. 

An alternative method to calculate r e ff is to directly integrate 
the He II Lya opacity from the CDDF. The only additional infor- 
mation needed is the distribution of line widths, provide d by t he 
Doppler parameter b. In this method, r c g is given by izuoll 19931) 



l + Z 



AHeII,Lyo: 

r N ii I, max 

x / d/VHi / dbf(Nm,b)W{Nni,b), (28) 

where W(Nn i, b) is the He II Lya equivalent width of an absorber 
with Doppler parameter b and f(Nni, b) is the joint distribution of 
A^h i and b. We assume that A^h i and b are uncorrected and that the 
distribution of b is a Dirac-delta function at b — 30 km s~ , a rep- 
resen tative approximation for Hi Lya forest systems fam et al] 
l200lh . In this method we do not subject the resulting optical depth 
to any extra normalization. 

The results of the FGPA and CDDF methods are shown in 
Figure [7] Both methods demonstrate that steep evolution of ThcH 
naturally leads to steep evolution in the observed r cr f . The addition 
of fluctuations further accelerates the evolution. The results for dif- 
ferent input parameters (a, Thi, CDDF) are largely the same (typ- 



ically < 20% at z ~ 3.4) in the FGPA method when normalized 
at z — 2. In contrast, the CDDF method depends sensitively on 
7/ t hi n oc Thi/ThcH (equation [6}, which can differ by a factor of 
a few between models. Thus, for a given CDDF, the He II optical 
depth places a joint constraint on Th i and a, subject to the uncer- 
tainties inherent in our cosmological radiative transfer model. 

For context, we also sho w measured r e g points in Figure [7] 
taken from the com pilation bv|Dixon & Furlanetto! J20 09) and the 
published results of IShull et all bold) andlWorseck et all l l201ll) . 
We re-binned the data from Worseck et al. I fcO 111) into Az = 0.12 
and 0.08 bins where appropriate and estimated lcr errors of the 
mean optical depth in each bin by combining the reported observa- 
tional errors with the observed scatter within the bin. Note that we 
have not averaged the recent data sets together because they appear 
to have systematic variations. For example, the data points from 
2 ~ 2.35-2.65 primarily represent measurements of the spectrum 
of the same quasar, HE 2347-4342, so the spread in their values be- 
tween the various groups does not represent true fluctuations in the 
optical depth but rather systematic differences in the data analysis. 
Those differences make it especially difficult to compare our mod- 
els to observations. Nevertheless, it is interesting that our models - 
which explicitly ignore He II reionization - match the evolution in 
the observed mean optical depth rather well. Additionally, the fluc- 
tuating background models appear to match the observations more 
closely than the uniform models, especially at z > 2.7 where the 
observed optical depth evolution is very steep. Our result demon- 
strates that the observed trend in and of itself does not require the 
He II fraction to evolve, although it also does not rule out such evo- 
lution. 

Unfortunately, our models do not explicitly describe how the 
integrated T cn - should vary at the same redshift along different 
lines of sight, even whe n averaged over large path lengths. Hydro- 
dynamic s imulations b y McOui nn et al.1 J2009I) and semi-analytic 
models by Furlanetto & Di xonl d2010h have described spatial varia- 
tions i n T e ff. Interestingly, the well-studied spectrum of HE 2347- 
4342 dReimers et all 1 19971 : iKriss et alj 1200 it IZheng et all |2004| : 
IShull et alj2004l20ld)" shows regions of high optical depth that ap- 
pear to require large swathes of He II at 2.7 < z < 2.9. We therefore 
emphasize that our models do not demand that He II reionization be 
over by z ~ 3; they instead demonstrate that, with respect to the 
evolution of the mean opacity, it is not required. 



5 DISCUSSION 

Our model for background fluctuations increases the average opac- 
ity of the IGM when the mean free path is comparable to the sep- 
aration between bright sources. This effect is primarily due to the 
skewness of f(T) towards lower V as the mean free path decreases 
(as in Figure [T}. While the absolute effect of our fluctuations pre- 
scription on the ionizing background may only be about a factor of 
a few, it predicts a steep increase in the ionizing background when 
the background transitions from being dominated by local sources 
to a smoother background with contributions from distant sources. 



5.1 Comparison to past theoretical work 

Figure [8] shows how our model c ompares to a pai r of r ecent 
ionizing background calculat ions by lHaardt & Madaul d2012l) and 
Faucher-Giguere et al. (2009). 



Faucher-Giguere et al. I feoosh used a single power-law CDDF 



with p — 1.4 and 7 = 1.5 that severely underestimates the num- 
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Figure 8. Uniform (solid black) and fluctuating (dashed black) He II ion- 
ization rate from this work com pared to the models fro m Haardt & Madaul 
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this is unlikely to be significant jPixon & Furlanettoll2012l) . The 
absorbers also show some clustering I Rudie et al 2012bl lah which 
will modulate the metagalactic radiation field (although likely only 
modestly). 

Other obvious sources of additional fluctuations in the ioniz- 
ing background - over and above those from the discrete sources 
- include radiative transfer effects (e.g. "shadows" behind op- 
tically thick regions as in iTittlev & Meiksinl l2007t) and c olli- 
sional ionization in superheated shocks dMuzahidetai]|201ll) . Of 
course, incomplete He II reionization may leave opaque patches 
of He I I that would introdu c e seve re fluctuations (McOui nn et alj 
120091; Furlanetto & Dixon 201 (3), possibly observed rec ently 
dZheng et alj |2004 IShull et alj l2010i: IWorseck et ail l201lh . We 
have explicitly ignored this possibility here so as to consider the 
evolution of the ionizing background in the absence of such effects. 

We also treat recombinations only approximately. We include 
recombination emission in our fluctuating model calculation in the 
same way as in the uniform model, by simply adding to the pre- 
existing quasars' emissivity. It therefore implicitly has the same 
source distribution, while in fact it will be more uniform than the 
point-like quasars because it is distributed throughout the IGM. On 
the other hand, recombination emission in low Thc ii regions will 
be weaker, and much of the emission from high Thc ii regions (i.e. 
near bright quasars) will not travel much beyond those quasar prox- 
imity regions before redshifting below ^Hcii ( < 30 Mpc ~ A m i„; 
ij |2.1.2t . so its effect on /(T) should be fairly minor. 



ber of low-density Lya forest absorbers compared to recent obser- 
vations (see the left panel of Figure [2} an d evolves more slowly 
than implied by Lya forest measurements dKim et alj|2002ah . Be- 
cause their CDDF severely underestimates the Hi opacity of the 
IGM from sub-LLS absorbers, they were forced to renormalize the 
quasar emissivity of ionizing photons at the hydrogen ionizing edge 
by a factor of 0.36 to match their measured Thi ~ 0.5 x 10~ 12 
s -1 llFaucher-Gi guere et alj2008al) . and thus their Thc ii is normal- 
ized somewhat lower as well. Their Fh c ii evolves at a similar rate 
to ou r fiducial uniform background model. 

lHaardt & Madaul d2012h used a CDDF that evolves more 
rapidly with redshift than our fiducial model and calculated a Th i 
that peaks at z ~ 2 and declines slowly towards higher redshift. 
They also used a different fitting form for the structure of IGM ab- 
sorbers. In their fit they more accurately approximated the average 
ionization rate within absorbers, which resulted in a more accu- 
rate fit to r\ at large H I column densities. However, as mentioned 
previously in § 12.1.11 those high Nh i systems do not contribute a 
substantial fraction of the opacity near the He II edge, and thus our 
approximation should not significantly affect our results. 

5.2 Fluctuating Model Caveats 

Other than the general simplifications necessary to invoke the cos- 
mological radiative transfer model, our parameterization of the 
fluctuations in the background is an ad hoc addition to a model 
designed for a medium with a uniform emissivity. In this section, 
we describe the primary uncertainties with such an approach. 

First of all, we may not have accurately captured the extent 
and character of the fluctuations. Spatial correlations in the ionizing 
background exist due to the large proximity regions of the primary 
sources (L» quasars). The massive hosts of these luminous quasars 
are clustered, which will increase the amplitude of the fluctua- 
tions. However, the proximity zones of the quasars are so large that 



6 CONCLUSION 

We have calculated the He II ionizing background using a cosmo- 
logical radiative transfer model that takes into account the lat- 
est constraints on quasar and IGM source properties. In our uni- 
form background model, which closely mimics previous work 
dFardaletal]|l998k iFaucher-Giguere et aljp2009l : lHaardt & Madaul 
|2012|) . we found that the He II ionization rate, Fh c ii, and the mean 
free path of He II ionizing photons should both evolve significantly 
during the time after Hell reionization (z ~ 2 — 3). However, at 
z ~ 3, the mean free path of He II ionizing photons is comparable 
to the average distance between the bright quasars that contribute 
most of the ionizing emissivity. While previous work investigated 
how this effect introduces fluctuations in th e ionizing backg round 
dFardal et alj|l99ct iMeiksin & White! 120031 : lFurlan"ettoll2009h , its 
implications for the mean ionizing background itself have not been 
studied in detail until now. 

We investigated for the first time how these fluctuations can 
affect the evolution of the mean background. We incorporated the 
distribution f(T) into our cosmological radiative transfer model 
by averaging the opacity to Hell ionizing photons over it. How- 
ever, that procedure still models the emission as diffuse sources 
rather than point-like quasars, so we supplemented it with a physi- 
cal model that accounts for the transparent zones around individual 
quasars when the sources become very rare. Including that model, 
our results showed that the fluctuating background introduces an- 
other source of opacity which causes the ionization rate to decrease 
by a factor of ~ 30% at z ~ 3.1 relative to the uniform background 
calculation. For z > 3.2 our transparent near-zone approximation 
dominates the results, so our model is no longer reliable - but we 
robustly predict very steep evolution in the ionizing background for 
a wide range of observed input parameters. 

As an example of the utility of our ionizing background 
model, we used the resulting ionization rate to estimate the evo- 
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lution of the Hell Lya effective optical depth, T e g. Rapid evolu- 
tion at z > 2.5, similar to that seen in observations, appears to be 
a natural consequence of a steeply evolving ionization rate. The 
addition of fluctuations improves our model's resemblance to the 
observed r e g evolution somewhat, though systematic uncertainties 
in the data analysis make a detailed comparison difficult. 

We note that our model does not incorporate 
He II reionization: that is, we assume that the He II fraction 
is very small throughout the IGM. We have therefore shown that 
reionization is not the only possible cause of a rapidly evolving 
ionizing background. Instead, the interaction between the (slowly) 
increasing emissivity and the (slowly) evolving IGM dumpiness 
can feed back on each other, strongly amplifying the evolution of 
the ionizing bac kground. Such evolu tion is naively predicted by 
simple models iMcOuinn et alj|201lh but is not observed in the 
hydrogen-ionizing background. 

Our result emphasizes the importance of understanding the 
IGM for interpreting measurements of the ionizing background and 
of reionization, includ ing that of both He II and H I . In the context 
of He II reionization, iDixon & Furlanetto! d2009l) argued that the 
rapidly increasing mean optical depth in the He II Lya line is con- 
sistent with ongoing Hell reionization at z>2.7. However, they 
prescribed a relatively slow evolution in the mean free path of ion- 
izing photons. 

On the other hand, a number of observations show substantial 
fluctuations in the mean optical depth, even when averaged over 
large scales. Our model does not address such large-scale fluc- 
tuations, because we have not incorporated any spatial informa- 
tion into the calculation. However, we do not expect th e conc lu- 
sions of iFurlanetto & Dixonl < |2010|) and IMcOuinn et all d2009h - 
that such fluctuations do not occur in a fully-ionized universe - to 
change substantially, because that is due to the shape of the proba- 
bility distribution of the ionizing background, and in particular the 
"floor" provided by distant sources. Thus it is very unlikely that 
even a small mean free path (~ 30 Mpc) for these photons can 
cause large opaque regions over several tens of Mpc, as obser ved 
dReimers et all 1 9971: IZheng et alj|2004l ; lshull et all2004ll2010h . 

This calculation may also have important implications for 
H I reionization, where an apparent rapid increase in the H I Lya 
optical depth has long b een attributed to the tail e nd of reionization 
{ Fan et al 1 I2002L l2006h . IFurlanetto & Mesinged (2009) previously 
showed that the overlap process of reionization (when ionized 
bubbles overlap to fill space) does not by itself cause a rapid 
increase in the ionizing background. We have shown that such an 
increase can be caused by "normal" post-reionization processes, 
through the interaction of a slowly increasing emissivity and 
slowly decreasing IGM clumping. Whether this occurs during 
Hi reionization cannot be said, because it depends sensitively 
on the evolution of that clumping (which is largely hidden due 
to the high opacity of the Lya forest beyond z ~ 6). However, 
this Hell analog indicates that a proper interpretation of data 
about H I reionization requires careful modelling (and ideally 
observations) of the IGM and not simply an understanding of the 
emitting sources. 

We thank K. Dixon, G. Rudie, and C. Steidel for helpful con- 
versations. This research was partially supported by the David and 
Lucile Packard Foundation and the Alfred R Sloan Foundation. 
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